Search for: All records

A three-dimensional, regional simulation is used to investigate ionospheric plasma density irregularities associated with Equatorial Spread F. This simulation is first driven with background electric fields derived from ISR observations. Next, the simulation is driven with electric fields taken from the WAM-IPE global model. The discrepancies between the two electric fields, particularly in the evening prereversal enhancement, produce disagreeing simulation results. The WAM-IPE electric fields are then studied through a simple sensitivity analysis of a field-line integrated electrodynamics model similar to the one used in WAM-IPE. This analysis suggests there is no simple tuning of ion composition or neutral winds that accurately reproduce ISR-observed electric fields on a day-to-day basis. Additionally, the persistency of the prereversal enhancement structure over time is studied and compared to measurements from the ICON satellite. These results suggest that WAM-IPE electric fields generally have a shorter and more variable correlation time than those measured by ICON. 

Low Earth orbit (LEO) radio occultation|radio occultations (RO) constellations can provide global electron density profiles (EDPs) to better specify and forecast the ionosphere‐thermosphere (I‐T) system. To inform future RO constellation design, this study uses comprehensive Observing System Simulation Experiments (OSSEs) to assess the ionospheric specification impact of assimilating synthetic EDPs into a coupled I‐T model. These OSSEs use 10 different sets of RO constellation configurations containing 6 or 12 LEO satellites with base orbit parameter combinations of 520 or 800 km altitude, and 24° or 72° inclination. The OSSEs are performed using the Ensemble Adjustment Kalman Filter implemented in the data assimilation (DA) Research Testbed and the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model (TIEGCM). A different I‐T model is used for the nature run, the Whole Atmosphere Model‐Ionosphere Plasmasphere Electrodynamics (WAM‐IPE), to simulate the period of interest is the St. Patrick's Day storm on March 13–18, 2015. Errors from models and EDP retrieval are realistically accounted for in this study through distinct I‐T models and by retrieving synthetic EDPs through an extension Abel inversion algorithm. OSSE assessment, using multiple metrics, finds that greater EDP spatial coverage leading to improved specification at altitudes 300 km and above, with the 520 km altitude constellations performing best due to yielding the highest observation counts. A potential performance limit is suggested with two 6‐satellite constellations. Lastly, close examination of Abel inversion error impacts highlights major EDP limitations at altitudes below 200 km and dayside equatorial regions with large horizontal gradients and low electron density magnitudes. 

This paper presents a study of the global medium‐scale (scales620 km) gravity wave (GW) activity (in terms of zonal wind variance) and its seasonal, local time, and longitudinal variations by employing the enhanced‐resolution (50 km) whole atmosphere model (WAMT254) and space‐based observations for geomagnetically quiet conditions. It is found that the GW hotspots produced by WAMT254 in the troposphere and stratosphere agree well with previously well‐studied orographic and nonorographic sources. In the ionosphere‐thermosphere (IT) region, GWs spread out forming latitudinal band‐like hotspots. During solstices, a primary maximum in GW activity is observed in WAMT254 and GOCE over winter mid‐high latitudes, likely associated with higher‐order waves with primary sources in polar night jet, fronts, and polar vortex. During all the seasons, the enhancement of GWs around the geomagnetic poles as observed by GOCE (at 250 km) is well captured by simulations. WAMT254 GWs in the IT region also show dependence on local time due to their interaction with migrating tides leading to diurnal and semidiurnal variations. The GWs are more likely to propagate up from the MLT region during westward/weakly eastward phase of thermospheric tides, signifying the dominance of eastward GW momentum flux in the MLT. Additionally, as a novel finding, a wavenumber‐4 signature in GW activity is predicted by WAMT254 between 6 and 12 local times in the tropics at 250 km, which propagates eastward with local time. This behavior is likely associated with the modulation of GWs by wave‐4 signal of nonmigrating tides in the lower thermospheric zonal winds. 

An expression of Rayleigh‐Taylor (R‐T) instability growth rate based on the field‐line integrated theory is newly established. This expression can be directly utilized in ionosphere models with magnetic flux tube structure based on Modified Apex Coordinates. The R‐T instability growth rates are calculated using the thermospheric and ionospheric conditions based on the coupled Whole Atmosphere Model and Ionosphere Plasmasphere Electrodynamic model (WAM‐IPE). The parameters used in this calculation include the field‐line integrated conductivities and currents, which consider the Quasi‐Dipole Coordinates and the modifications to the equations of electrodynamics. Detailed description of the new formulas and comprehensive analyses of diurnal, longitudinal, and seasonal variations of the R‐T instability growth rate are carried out. The dependencies of growth rates on pre‐reversal enhancement (PRE) vertical drifts and solar activity are also examined. The results show that pronounced R‐T growth rates are captured between 18 and 22 local time (LT) when strong PRE occurs in the equatorial ionosphere. The simulated R‐T growth rate increases with increasing solar activity levels and demonstrates strong correlations with the angle between the sunset terminator and the geomagnetic field line. These results are consistent with plasma irregularity occurrence rates shown in various satellite observations, suggesting that the newly developed R‐T growth rate calculation can effectively capture the probability of irregularities by considering the changes along magnetic flux‐tubes in the ionosphere. Since the WAM‐IPE is running in operation at National Oceanic and Atmospheric Administration Space Weather Prediction Center, the new calculations can be potentially implemented in the near future to provide forecasted information of the R‐T growth rate. 

Key Points Validation of ionospheric total electron content (TEC) by the state‐of‐the‐art ionospheric models hosted by NASA Community Coordinated Modeling Center, National Oceanic and Atmospheric Administration Space Weather Prediction Center, and NASA Jet Propulsion Laboratory (JPL) Multiple metrics and skill scores are used to assess the performance of ionospheric models in capturing storm time TEC anomaly GLObal Total Electron Content and JPL Global Ionospheric Map perform best, and physics‐based models perform better than the empirical model in capturing storm TEC variations 

Abstract This study presents a data‐driven approach to quantify uncertainties in the ionosphere‐thermosphere (IT) system due to varying solar wind parameters (drivers) during quiet conditions (Kp < 4) and fixed solar radiation and lower atmospheric conditions representative of 16 March 2013. Ensemble simulations of the coupled Whole Atmosphere Model with Ionosphere Plasmasphere Electrodynamics (WAM‐IPE) driven by synthetic solar wind drivers generated through a multi‐channel variational autoencoder (MCVAE) model are obtained. Applying the polynomial chaos expansion (PCE) technique, it is possible to estimate the means and variances of the QoIs as well as the sensitivities of the QoIs with regard to the drivers. Our results highlight unique features of the IT system's uncertainty: (a) the uncertainty of the IT system is larger during nighttime; (b) the spatial distributions of the uncertainty for electron density and zonal drift at fixed local times present 4 peaks in the evening sector, which are associated with the low‐density regions of longitude structure of electron density; (c) the uncertainty of the equatorial electron density is highly correlated with the uncertainty of the zonal drift, especially in the evening sector, while it is weakly correlated with the vertical drift. A variance‐based global sensitivity analysis suggests that the IMF Bz plays a dominant role in the uncertainty of electron density. A further discussion shows that the uncertainty of the IT system is determined by the magnitudes and universal time variations of solar wind drivers. Its temporal and spatial distribution can be modulated by the average state of the IT system. 

Abstract The Iterative Driver Estimation and Assimilation (IDEA) data assimilation technique was used with the Whole Atmosphere Model (WAM) to improve neutral density specification in the upper thermosphere. Two different neutral density data sources were examined to enhance the capability of simulating the global thermospheric state. The first were accelerometer estimates of neutral density from the Challenging Mini‐Satellite Payload (CHAMP) satellite. The second were neutral density estimates from the Global Ultraviolet Imager (GUVI) limb‐scan airglow observations aboard the Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite. Due to the intensity of the November 2003 storm, two changes were necessary in WAM. The first was allowing the Kp geomagnetic index to exceed 9 and the second was changing the relationship between Kp and the solar wind parameters used to drive the model. With these changes, results show that IDEA effectively captures the thermospheric neutral density at the CHAMP satellite altitude and follows the time‐dependence through the November 2003 storm period. Furthermore, a cross‐comparison was conducted with the GUVI dayside limb scan measurements. GUVI neutral densities within 270–320 km show the closest agreement with WAM when CHAMP data was assimilated by IDEA. We speculate on the potential for observations from GUVI at 300 km to be used as a data source in the IDEA‐WAM simulations. These simulations demonstrate the utility of the IDEA data assimilation technique with physical models and that using either accelerometer observations or ultraviolet airglow limb measurement during extreme storm periods could be used. 

Abstract Spread‐F (SF) is a feature that can be visually observed on ionograms when the ionosonde signals are significantly impacted by plasma irregularities in the ionosphere. Depending on the scale of the plasma irregularities, radio waves of different frequencies are impacted differently when the signals pass through the ionosphere. An automated method for detecting SF in ionograms is presented in this study. Through detecting the existence of SF in ionograms, we can help identify instances of plasma irregularities that are potentially affecting the high‐frequency radio‐wave systems. The ionogram images from Jicamarca observatory in Peru, during the years 2008–2019, are used in this study. Three machine learning approaches have been carried out: supervised learning using Support Vector Machines, and two neural network‐based learning methods: autoencoder and transfer learning. Of these three methods, the transfer learning approach, which uses convolutional neural network architectures, demonstrates the best performance. The best existing architecture that is suitable for this problem appears to be the ResNet50. With respect to the training epoch number, the ResNet50 showed the greatest change in the metric values for the key metrics that we were tracking. Furthermore, on a test set of 2050 ionograms, the model based on the ResNet50 architecture provides an accuracy of 89%, recall of 87%, precision of 95%, as well as Area Under the Curve of 96%. The work also provides a labeled data set of around 28,000 ionograms, which is extremely useful for the community for future machine learning studies.